------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  /dataverse files/Re-Assessing Governmental Corruption of Religion-2026-logfile.log
  log type:  text
 opened on:  29 Jan 2026, 09:54:06

. * Re-Assessing Governmental Corruption of Religion: Where Political Theory Meets Empirical Evidence *

. * Replication DO File *

. 
. *********************************

. * All Christian Majority States *

. *********************************

. asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class pol
> itics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, nested dec(
> 3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) reset

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85630.315  
Iteration 2:  Log likelihood = -85372.594  
Iteration 3:  Log likelihood = -85371.867  
Iteration 4:  Log likelihood = -85371.867  

Refining starting values:

Grid node 0:  Log likelihood = -82754.802

Fitting full model:

Iteration 0:  Log likelihood = -82754.802  (not concave)
Iteration 1:  Log likelihood = -82748.023  (not concave)
Iteration 2:  Log likelihood = -82741.652  (backed up)
Iteration 3:  Log likelihood = -82669.758  
Iteration 4:  Log likelihood = -82569.519  
Iteration 5:  Log likelihood = -82469.532  
Iteration 6:  Log likelihood =  -82441.29  
Iteration 7:  Log likelihood = -82430.806  
Iteration 8:  Log likelihood =  -82430.75  
Iteration 9:  Log likelihood =  -82430.75  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   11989.14
Log likelihood = -82430.75                      Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
     lz06funding_l1 |   .0281889   .0139422     2.02   0.043     .0008626    .0555151
             mxx_l1 |   .0058173   .0030766     1.89   0.059    -.0002127    .0118474
      church_attend |    .288619   .0037609    76.74   0.000     .2812479    .2959902
          religious |   .8777018   .0198287    44.26   0.000     .8388383    .9165654
            atheist |  -.6358739   .0699162    -9.09   0.000    -.7729071   -.4988407
     satisfied_life |   .0311928    .003349     9.31   0.000     .0246289    .0377567
              edu08 |  -.0913121   .0035866   -25.46   0.000    -.0983417   -.0842825
       social_class |  -.0132357   .0081307    -1.63   0.104    -.0291715    .0027001
       politics_l2r |   .0471544   .0031741    14.86   0.000     .0409333    .0533754
          zpcgdplog |   .1628553   .0771961     2.11   0.035     .0115538    .3141569
           zpolity2 |   .0031426   .0049527     0.63   0.526    -.0065644    .0128497
           zdurable |  -.0078576   .0026374    -2.98   0.003    -.0130268   -.0026884
          ezmajpctx |  -.0178103   .0046691    -3.81   0.000    -.0269615   -.0086592
            zpoplog |  -.2454276   .0957104    -2.56   0.010    -.4330165   -.0578387
    fraction_tc_min |  -3.083792    .417558    -7.39   0.000    -3.902191   -2.265394
                    |
             region |
       Former USSR  |  -.6844166   .3845175    -1.78   0.075    -1.438057     .069224
              Asia  |   1.391645   .9954373     1.40   0.162     -.559376    3.342666
Sub Saharan Africa  |   .9439547   .4651709     2.03   0.042     .0322366    1.855673
     Latin America  |   .1249988   .3930263     0.32   0.750    -.6453186    .8953163
--------------------+----------------------------------------------------------------
              /cut1 |  -6.940299   1.939921                     -10.74247   -3.138124
              /cut2 |  -4.664661   1.939713                      -8.46643   -.8628929
              /cut3 |  -2.622667   1.939592                     -6.424198    1.178864
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .8665677   .2034077                      .5470181    1.372787
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 5882.23       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class polit
> ics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, nested dec(3)
>  tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85578.777  
Iteration 2:  Log likelihood = -85320.969  
Iteration 3:  Log likelihood = -85320.236  
Iteration 4:  Log likelihood = -85320.236  

Refining starting values:

Grid node 0:  Log likelihood = -82751.627

Fitting full model:

Iteration 0:  Log likelihood = -82751.627  (not concave)
Iteration 1:  Log likelihood = -82744.797  (not concave)
Iteration 2:  Log likelihood = -82741.632  (backed up)
Iteration 3:  Log likelihood = -82667.516  
Iteration 4:  Log likelihood = -82565.337  
Iteration 5:  Log likelihood = -82462.116  
Iteration 6:  Log likelihood = -82432.373  
Iteration 7:  Log likelihood = -82421.263  
Iteration 8:  Log likelihood = -82421.188  
Iteration 9:  Log likelihood = -82421.188  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   12005.87
Log likelihood = -82421.188                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
       vfundingx_l1 |  -.2543299   .0540828    -4.70   0.000    -.3603302   -.1483297
             mxx_l1 |   .0068141   .0030977     2.20   0.028     .0007428    .0128854
      church_attend |   .2885299    .003761    76.72   0.000     .2811584    .2959014
          religious |   .8761553   .0198312    44.18   0.000     .8372868    .9150237
            atheist |  -.6417827   .0699076    -9.18   0.000    -.7787991   -.5047664
     satisfied_life |   .0312774   .0033489     9.34   0.000     .0247138    .0378411
              edu08 |  -.0912809   .0035869   -25.45   0.000    -.0983111   -.0842507
       social_class |  -.0138862   .0081324    -1.71   0.088    -.0298254    .0020529
       politics_l2r |   .0470009   .0031742    14.81   0.000     .0407796    .0532223
          zpcgdplog |   .1616952   .0767292     2.11   0.035     .0113087    .3120817
           zpolity2 |   .0036558   .0049453     0.74   0.460    -.0060368    .0133484
           zdurable |  -.0068859   .0026172    -2.63   0.009    -.0120154   -.0017563
          ezmajpctx |  -.0180203   .0045659    -3.95   0.000    -.0269693   -.0090712
            zpoplog |  -.2302189      .0995    -2.31   0.021    -.4252353   -.0352026
    fraction_tc_min |    -3.2902   .4121777    -7.98   0.000    -4.098054   -2.482347
                    |
             region |
       Former USSR  |  -.4243227    .392082    -1.08   0.279    -1.192789     .344144
              Asia  |   1.463924   1.040657     1.41   0.160    -.5757258    3.503574
Sub Saharan Africa  |   .6737275   .4890047     1.38   0.168    -.2847042    1.632159
     Latin America  |   .3066156   .4098282     0.75   0.454     -.496633    1.109864
--------------------+----------------------------------------------------------------
              /cut1 |  -7.533581   2.001203                     -11.45587   -3.611296
              /cut2 |  -5.257719   2.001004                     -9.179615   -1.335823
              /cut3 |  -3.215087   2.000873                     -7.136727    .7065522
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .9528289   .2207482                      .6050781    1.500439
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 5798.10       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_clas
> s politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, nested
>  dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85642.458  
Iteration 2:  Log likelihood = -85385.817  
Iteration 3:  Log likelihood = -85385.081  
Iteration 4:  Log likelihood = -85385.081  

Refining starting values:

Grid node 0:  Log likelihood =  -82748.28

Fitting full model:

Iteration 0:  Log likelihood =  -82748.28  (not concave)
Iteration 1:  Log likelihood = -82741.473  (backed up)
Iteration 2:  Log likelihood = -82736.726  (backed up)
Iteration 3:  Log likelihood = -82664.675  
Iteration 4:  Log likelihood = -82565.153  
Iteration 5:  Log likelihood = -82465.497  
Iteration 6:  Log likelihood = -82437.267  
Iteration 7:  Log likelihood = -82426.798  
Iteration 8:  Log likelihood = -82426.741  
Iteration 9:  Log likelihood = -82426.741  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   11995.30
Log likelihood = -82426.741                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |  -.0927677   .0266623    -3.48   0.001    -.1450249   -.0405104
             mxx_l1 |   .0045935   .0030939     1.48   0.138    -.0014705    .0106575
      church_attend |    .288736    .003761    76.77   0.000     .2813646    .2961073
          religious |   .8769958   .0198294    44.23   0.000     .8381308    .9158607
            atheist |   -.645832   .0699536    -9.23   0.000    -.7829386   -.5087254
     satisfied_life |   .0311186   .0033487     9.29   0.000     .0245552     .037682
              edu08 |  -.0915676   .0035877   -25.52   0.000    -.0985994   -.0845359
       social_class |  -.0137969   .0081319    -1.70   0.090     -.029735    .0021413
       politics_l2r |   .0470076   .0031743    14.81   0.000     .0407861    .0532291
          zpcgdplog |   .2238082   .0761193     2.94   0.003     .0746172    .3729993
           zpolity2 |   .0024269   .0049543     0.49   0.624    -.0072833     .012137
           zdurable |  -.0072119   .0025872    -2.79   0.005    -.0122827   -.0021412
          ezmajpctx |  -.0157509   .0046241    -3.41   0.001     -.024814   -.0066878
            zpoplog |  -.2814766   .0985131    -2.86   0.004    -.4745587   -.0883946
    fraction_tc_min |  -3.194759   .4174929    -7.65   0.000     -4.01303   -2.376488
                    |
             region |
       Former USSR  |  -.4581258    .381864    -1.20   0.230    -1.206566     .290314
              Asia  |   1.649603   1.002732     1.65   0.100    -.3157154    3.614922
Sub Saharan Africa  |   1.207318   .4690428     2.57   0.010      .288011    2.126625
     Latin America  |   .2129496   .3949811     0.54   0.590    -.5611991    .9870983
--------------------+----------------------------------------------------------------
              /cut1 |  -6.939627   1.981942                     -10.82416   -3.055092
              /cut2 |  -4.663633   1.981733                     -8.547758   -.7795074
              /cut3 |  -2.621301   1.981599                     -6.505164    1.262562
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .8819603    .207134                      .5565945    1.397524
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 5916.68       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class pol
> itics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, nested dec(
> 3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood = -18949.878  
Iteration 2:  Log likelihood = -18933.967  
Iteration 3:  Log likelihood = -18933.944  
Iteration 4:  Log likelihood = -18933.944  

Refining starting values:

Grid node 0:  Log likelihood = -18157.202

Fitting full model:

Iteration 0:  Log likelihood = -18157.202  
Iteration 1:  Log likelihood = -18136.124  
Iteration 2:  Log likelihood = -18092.148  
Iteration 3:  Log likelihood = -18074.328  
Iteration 4:  Log likelihood = -18074.038  
Iteration 5:  Log likelihood = -18074.033  
Iteration 6:  Log likelihood = -18074.033  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1616.00
Log likelihood = -18074.033                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
     lz06funding_l1 |  -.1001633   .0258016    -3.88   0.000    -.1507334   -.0495931
             mxx_l1 |  -.0171478   .0081811    -2.10   0.036    -.0331824   -.0011131
      church_attend |   .2136186   .0075845    28.17   0.000     .1987532     .228484
          religious |   .6381599   .0441664    14.45   0.000     .5515953    .7247244
            atheist |  -.2201592   .1548197    -1.42   0.155    -.5236003    .0832818
     satisfied_life |   .0439243   .0072548     6.05   0.000     .0297051    .0581435
              edu08 |  -.0725832     .00789    -9.20   0.000    -.0880473    -.057119
       social_class |   .0010525   .0163898     0.06   0.949    -.0310708    .0331759
       politics_l2r |   .0494329   .0066979     7.38   0.000     .0363054    .0625605
          zpcgdplog |  -.1241533   .1433751    -0.87   0.387    -.4051633    .1568567
           zpolity2 |  -.0377301   .0120816    -3.12   0.002    -.0614096   -.0140506
           zdurable |   .0032531   .0043667     0.74   0.456    -.0053055    .0118117
          ezmajpctx |  -.0020879   .0104645    -0.20   0.842    -.0225979    .0184222
            zpoplog |  -.0238505   .1188376    -0.20   0.841     -.256768    .2090669
    fraction_tc_min |  -2.915481    1.09117    -2.67   0.008    -5.054136   -.7768263
                    |
             region |
       Former USSR  |   .1814228    .516278     0.35   0.725    -.8304635    1.193309
              Asia  |   .7023653    1.15021     0.61   0.541    -1.552004    2.956735
Sub Saharan Africa  |   .7315643   .6384826     1.15   0.252    -.5198385    1.982967
     Latin America  |  -.4540791   .4841978    -0.94   0.348    -1.403089    .4949311
--------------------+----------------------------------------------------------------
              /cut1 |  -4.457378   3.050337                     -10.43593    1.521173
              /cut2 |  -2.838361   3.049647                     -8.815559    3.138837
              /cut3 |  -1.108958   3.049097                     -7.085077    4.867162
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   1.049933   .2774822                      .6254623     1.76247
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 1719.82       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class polit
> ics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, nested dec(3)
>  tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood =  -18951.51  
Iteration 2:  Log likelihood = -18935.781  
Iteration 3:  Log likelihood = -18935.759  
Iteration 4:  Log likelihood = -18935.759  

Refining starting values:

Grid node 0:  Log likelihood = -18161.658

Fitting full model:

Iteration 0:  Log likelihood = -18161.658  
Iteration 1:  Log likelihood =  -18136.07  
Iteration 2:  Log likelihood = -18084.975  
Iteration 3:  Log likelihood = -18062.489  
Iteration 4:  Log likelihood = -18061.514  
Iteration 5:  Log likelihood = -18061.468  
Iteration 6:  Log likelihood = -18061.468  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1631.36
Log likelihood = -18061.468                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
       vfundingx_l1 |  -.4518231   .0758368    -5.96   0.000    -.6004605   -.3031857
             mxx_l1 |  -.0096081   .0084999    -1.13   0.258    -.0262677    .0070514
      church_attend |    .213148   .0075838    28.11   0.000      .198284     .228012
          religious |   .6377505   .0441784    14.44   0.000     .5511624    .7243385
            atheist |  -.2279811   .1548875    -1.47   0.141    -.5315551    .0755929
     satisfied_life |   .0439635   .0072536     6.06   0.000     .0297468    .0581802
              edu08 |  -.0727266    .007896    -9.21   0.000    -.0882024   -.0572508
       social_class |   .0004156   .0163952     0.03   0.980    -.0317184    .0325495
       politics_l2r |   .0489411    .006701     7.30   0.000     .0358073    .0620748
          zpcgdplog |  -.3808488   .1538891    -2.47   0.013    -.6824658   -.0792318
           zpolity2 |  -.0392709    .012115    -3.24   0.001    -.0630159   -.0155259
           zdurable |   .0081499   .0047628     1.71   0.087     -.001185    .0174849
          ezmajpctx |  -.0114586   .0109129    -1.05   0.294    -.0328474    .0099302
            zpoplog |   .0142215   .1265469     0.11   0.911    -.2338059    .2622488
    fraction_tc_min |  -4.418305   1.151241    -3.84   0.000    -6.674696   -2.161913
                    |
             region |
       Former USSR  |   .1043013   .5394173     0.19   0.847    -.9529371     1.16154
              Asia  |   .2307314   1.240303     0.19   0.852    -2.200219    2.661682
Sub Saharan Africa  |   .0561877   .6756982     0.08   0.934    -1.268156    1.380532
     Latin America  |  -.2251377   .5189032    -0.43   0.664    -1.242169    .7918939
--------------------+----------------------------------------------------------------
              /cut1 |  -8.444036   3.246118                     -14.80631   -2.081762
              /cut2 |  -6.822377   3.245432                     -13.18331   -.4614475
              /cut3 |  -5.090229   3.244721                     -11.44976    1.269308
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|    1.24168   .3238958                      .7446828     2.07037
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 1748.58       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_clas
> s politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, nested
>  dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood = -18943.101  
Iteration 2:  Log likelihood = -18927.443  
Iteration 3:  Log likelihood =  -18927.42  
Iteration 4:  Log likelihood =  -18927.42  

Refining starting values:

Grid node 0:  Log likelihood = -18159.397

Fitting full model:

Iteration 0:  Log likelihood = -18159.397  
Iteration 1:  Log likelihood = -18147.577  
Iteration 2:  Log likelihood =   -18117.6  
Iteration 3:  Log likelihood = -18087.999  
Iteration 4:  Log likelihood = -18079.613  
Iteration 5:  Log likelihood = -18076.489  
Iteration 6:  Log likelihood =  -18076.48  
Iteration 7:  Log likelihood =  -18076.48  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1614.85
Log likelihood = -18076.48                      Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |   .2752741   .0849802     3.24   0.001     .1087159    .4418323
             mxx_l1 |  -.0214377    .007962    -2.69   0.007    -.0370429   -.0058324
      church_attend |   .2147902   .0075806    28.33   0.000     .1999324     .229648
          religious |   .6394209   .0441713    14.48   0.000     .5528466    .7259951
            atheist |  -.2215843   .1548086    -1.43   0.152    -.5250036    .0818351
     satisfied_life |   .0450795   .0072497     6.22   0.000     .0308704    .0592886
              edu08 |  -.0719142   .0078849    -9.12   0.000    -.0873684     -.05646
       social_class |   .0015925   .0163861     0.10   0.923    -.0305237    .0337086
       politics_l2r |   .0498847   .0066963     7.45   0.000     .0367603    .0630091
          zpcgdplog |  -.2110344   .1367045    -1.54   0.123    -.4789703    .0569016
           zpolity2 |  -.0406334   .0120273    -3.38   0.001    -.0642065   -.0170603
           zdurable |  -.0016355    .003799    -0.43   0.667    -.0090813    .0058103
          ezmajpctx |   .0058819   .0100988     0.58   0.560    -.0139113    .0256752
            zpoplog |   .0016194   .1113965     0.01   0.988    -.2167137    .2199525
    fraction_tc_min |  -1.592802   1.026256    -1.55   0.121    -3.604228    .4186231
                    |
             region |
       Former USSR  |   -.330756   .4954609    -0.67   0.504    -1.301841    .6403294
              Asia  |   .1734151   1.091148     0.16   0.874    -1.965197    2.312027
Sub Saharan Africa  |   .1568907   .6162192     0.25   0.799    -1.050877    1.364658
     Latin America  |  -.6232918   .4624896    -1.35   0.178    -1.529755    .2831711
--------------------+----------------------------------------------------------------
              /cut1 |  -3.500744   2.794174                     -8.977225    1.975736
              /cut2 |  -1.881682   2.793629                     -7.357095     3.59373
              /cut3 |   -.153785   2.793213                     -5.628382    5.320812
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .9312984   .2321379                      .5713683    1.517964
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 1701.88       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. **********************************************

. * All Christian Majority States, Odds Ratios *

. **********************************************

. asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class pol
> itics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, or nested d
> ec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) reset
(File Myfile.doc already exists, option append was assumed)

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85630.315  
Iteration 2:  Log likelihood = -85372.594  
Iteration 3:  Log likelihood = -85371.867  
Iteration 4:  Log likelihood = -85371.867  

Refining starting values:

Grid node 0:  Log likelihood = -82754.802

Fitting full model:

Iteration 0:  Log likelihood = -82754.802  (not concave)
Iteration 1:  Log likelihood = -82748.023  (not concave)
Iteration 2:  Log likelihood = -82741.652  (backed up)
Iteration 3:  Log likelihood = -82669.758  
Iteration 4:  Log likelihood = -82569.519  
Iteration 5:  Log likelihood = -82469.532  
Iteration 6:  Log likelihood =  -82441.29  
Iteration 7:  Log likelihood = -82430.806  
Iteration 8:  Log likelihood =  -82430.75  
Iteration 9:  Log likelihood =  -82430.75  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   11989.14
Log likelihood = -82430.75                      Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Odds ratio   Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
     lz06funding_l1 |    1.02859   .0143408     2.02   0.043     1.000863    1.057085
             mxx_l1 |   1.005834   .0030946     1.89   0.059     .9997873    1.011918
      church_attend |   1.334583   .0050192    76.74   0.000     1.324782    1.344457
          religious |   2.405365   .0476953    44.26   0.000     2.313678    2.500687
            atheist |   .5294726   .0370187    -9.09   0.000      .461669    .6072342
     satisfied_life |   1.031684   .0034551     9.31   0.000     1.024935    1.038479
              edu08 |   .9127328   .0032736   -25.46   0.000     .9063391    .9191716
       social_class |   .9868515   .0080238    -1.63   0.104     .9712499    1.002704
       politics_l2r |   1.048284   .0033273    14.86   0.000     1.041783    1.054826
          zpcgdplog |   1.176866   .0908495     2.11   0.035     1.011621    1.369105
           zpolity2 |   1.003148   .0049683     0.63   0.526     .9934571    1.012933
           zdurable |   .9921732   .0026168    -2.98   0.003     .9870577    .9973153
          ezmajpctx |   .9823473   .0045866    -3.81   0.000     .9733987    .9913782
            zpoplog |   .7823699   .0748809    -2.56   0.010     .6485498    .9438021
    fraction_tc_min |   .0457853    .019118    -7.39   0.000     .0201976    .1037892
                    |
             region |
       Former USSR  |   .5043844   .1939447    -1.78   0.075     .2373885    1.071676
              Asia  |   4.021461   4.003112     1.40   0.162     .5715656    28.29447
Sub Saharan Africa  |   2.570126   1.195548     2.03   0.042     1.032762    6.396001
     Latin America  |   1.133147   .4453566     0.32   0.750     .5244954     2.44811
--------------------+----------------------------------------------------------------
              /cut1 |  -6.940299   1.939921                     -10.74247   -3.138124
              /cut2 |  -4.664661   1.939713                      -8.46643   -.8628929
              /cut3 |  -2.622667   1.939592                     -6.424198    1.178864
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .8665677   .2034077                      .5470181    1.372787
-------------------------------------------------------------------------------------
Note: Estimates are transformed only in the first equation to odds ratios.
LR test vs. ologit model: chibar2(01) = 5882.23       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class polit
> ics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, or nested dec
> (3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85578.777  
Iteration 2:  Log likelihood = -85320.969  
Iteration 3:  Log likelihood = -85320.236  
Iteration 4:  Log likelihood = -85320.236  

Refining starting values:

Grid node 0:  Log likelihood = -82751.627

Fitting full model:

Iteration 0:  Log likelihood = -82751.627  (not concave)
Iteration 1:  Log likelihood = -82744.797  (not concave)
Iteration 2:  Log likelihood = -82741.632  (backed up)
Iteration 3:  Log likelihood = -82667.516  
Iteration 4:  Log likelihood = -82565.337  
Iteration 5:  Log likelihood = -82462.116  
Iteration 6:  Log likelihood = -82432.373  
Iteration 7:  Log likelihood = -82421.263  
Iteration 8:  Log likelihood = -82421.188  
Iteration 9:  Log likelihood = -82421.188  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   12005.87
Log likelihood = -82421.188                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Odds ratio   Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
       vfundingx_l1 |   .7754359   .0419377    -4.70   0.000      .697446    .8621469
             mxx_l1 |   1.006837   .0031188     2.20   0.028     1.000743    1.012969
      church_attend |   1.334464    .005019    76.72   0.000     1.324663    1.344338
          religious |   2.401648   .0476276    44.18   0.000     2.310091    2.496835
            atheist |   .5263532   .0367961    -9.18   0.000     .4589569    .6036466
     satisfied_life |   1.031772   .0034553     9.34   0.000     1.025022    1.038566
              edu08 |   .9127613    .003274   -25.45   0.000     .9063669    .9192008
       social_class |   .9862098   .0080202    -1.71   0.088      .970615    1.002055
       politics_l2r |   1.048123    .003327    14.81   0.000     1.041623    1.054664
          zpcgdplog |   1.175502   .0901953     2.11   0.035     1.011373    1.366266
           zpolity2 |   1.003663   .0049634     0.74   0.460     .9939814    1.013438
           zdurable |   .9931378   .0025992    -2.63   0.009     .9880565    .9982452
          ezmajpctx |   .9821411   .0044844    -3.95   0.000     .9733911    .9909698
            zpoplog |   .7943597   .0790388    -2.31   0.021      .653616    .9654098
    fraction_tc_min |   .0372464   .0153521    -7.98   0.000      .016605    .0835469
                    |
             region |
       Former USSR  |   .6542128   .2565051    -1.08   0.279     .3033739    1.410782
              Asia  |   4.322889   4.498643     1.41   0.160     .5622966      33.234
Sub Saharan Africa  |   1.961535   .9592001     1.38   0.168     .7522368    5.114907
     Latin America  |   1.358819   .5568822     0.75   0.454     .6085763    3.033946
--------------------+----------------------------------------------------------------
              /cut1 |  -7.533581   2.001203                     -11.45587   -3.611296
              /cut2 |  -5.257719   2.001004                     -9.179615   -1.335823
              /cut3 |  -3.215087   2.000873                     -7.136727    .7065522
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .9528289   .2207482                      .6050781    1.500439
-------------------------------------------------------------------------------------
Note: Estimates are transformed only in the first equation to odds ratios.
LR test vs. ologit model: chibar2(01) = 5798.10       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_clas
> s politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:, or nes
> ted dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85642.458  
Iteration 2:  Log likelihood = -85385.817  
Iteration 3:  Log likelihood = -85385.081  
Iteration 4:  Log likelihood = -85385.081  

Refining starting values:

Grid node 0:  Log likelihood =  -82748.28

Fitting full model:

Iteration 0:  Log likelihood =  -82748.28  (not concave)
Iteration 1:  Log likelihood = -82741.473  (backed up)
Iteration 2:  Log likelihood = -82736.726  (backed up)
Iteration 3:  Log likelihood = -82664.675  
Iteration 4:  Log likelihood = -82565.153  
Iteration 5:  Log likelihood = -82465.497  
Iteration 6:  Log likelihood = -82437.267  
Iteration 7:  Log likelihood = -82426.798  
Iteration 8:  Log likelihood = -82426.741  
Iteration 9:  Log likelihood = -82426.741  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   11995.30
Log likelihood = -82426.741                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Odds ratio   Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |   .9114052   .0243002    -3.48   0.001     .8650008    .9602992
             mxx_l1 |   1.004604   .0031082     1.48   0.138     .9985305    1.010715
      church_attend |   1.334739   .0050199    76.77   0.000     1.324937    1.344614
          religious |   2.403668   .0476633    44.23   0.000     2.312041    2.498925
            atheist |   .5242262   .0366715    -9.23   0.000     .4570609    .6012615
     satisfied_life |   1.031608   .0034546     9.29   0.000     1.024859    1.038401
              edu08 |   .9124996   .0032738   -25.52   0.000     .9061056    .9189387
       social_class |   .9862979   .0080204    -1.70   0.090     .9707027    1.002144
       politics_l2r |    1.04813   .0033271    14.81   0.000     1.041629    1.054671
          zpcgdplog |   1.250831   .0952124     2.94   0.003     1.077472    1.452083
           zpolity2 |    1.00243   .0049663     0.49   0.624     .9927431    1.012211
           zdurable |    .992814   .0025686    -2.79   0.005     .9877924    .9978611
          ezmajpctx |   .9843725   .0045519    -3.41   0.001     .9754913    .9933345
            zpoplog |   .7546686   .0743447    -2.86   0.004     .6221596    .9153996
    fraction_tc_min |   .0409764   .0171073    -7.65   0.000     .0180785    .0928762
                    |
             region |
       Former USSR  |   .6324679   .2415168    -1.20   0.230     .2992232    1.336847
              Asia  |   5.204915   5.219135     1.65   0.100      .729267    37.14845
Sub Saharan Africa  |   3.344503   1.568715     2.57   0.010     1.333772    8.386515
     Latin America  |   1.237322   .4887189     0.54   0.590     .5705246    2.683437
--------------------+----------------------------------------------------------------
              /cut1 |  -6.939627   1.981942                     -10.82416   -3.055092
              /cut2 |  -4.663633   1.981733                     -8.547758   -.7795074
              /cut3 |  -2.621301   1.981599                     -6.505164    1.262562
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .8819603    .207134                      .5565945    1.397524
-------------------------------------------------------------------------------------
Note: Estimates are transformed only in the first equation to odds ratios.
LR test vs. ologit model: chibar2(01) = 5916.68       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class pol
> itics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, or nested d
> ec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood = -18949.878  
Iteration 2:  Log likelihood = -18933.967  
Iteration 3:  Log likelihood = -18933.944  
Iteration 4:  Log likelihood = -18933.944  

Refining starting values:

Grid node 0:  Log likelihood = -18157.202

Fitting full model:

Iteration 0:  Log likelihood = -18157.202  
Iteration 1:  Log likelihood = -18136.124  
Iteration 2:  Log likelihood = -18092.148  
Iteration 3:  Log likelihood = -18074.328  
Iteration 4:  Log likelihood = -18074.038  
Iteration 5:  Log likelihood = -18074.033  
Iteration 6:  Log likelihood = -18074.033  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1616.00
Log likelihood = -18074.033                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Odds ratio   Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
     lz06funding_l1 |   .9046897   .0233424    -3.88   0.000      .860077    .9516165
             mxx_l1 |   .9829984    .008042    -2.10   0.036     .9673621    .9988875
      church_attend |    1.23815   .0093908    28.17   0.000     1.219881    1.256693
          religious |   1.892994   .0836067    14.45   0.000      1.73602    2.064162
            atheist |    .802391   .1242259    -1.42   0.155     .5923839    1.086848
     satisfied_life |   1.044903   .0075806     6.05   0.000     1.030151    1.059867
              edu08 |   .9299884   .0073376    -9.20   0.000     .9157176    .9444817
       social_class |   1.001053    .016407     0.06   0.949     .9694069    1.033732
       politics_l2r |   1.050675   .0070373     7.38   0.000     1.036972    1.064559
          zpcgdplog |   .8832444   .1266353    -0.87   0.387     .6668679    1.169828
           zpolity2 |   .9629728   .0116343    -3.12   0.002     .9404379    .9860476
           zdurable |   1.003258    .004381     0.74   0.456     .9947085    1.011882
          ezmajpctx |   .9979143   .0104427    -0.20   0.842     .9776555    1.018593
            zpoplog |   .9764316   .1160368    -0.20   0.841     .7735476    1.232527
    fraction_tc_min |    .054178   .0591174    -2.67   0.008     .0063829    .4598632
                    |
             region |
       Former USSR  |   1.198922   .6189771     0.35   0.725     .4358472    3.297977
              Asia  |   2.018522   2.321723     0.61   0.541     .2118231    19.23506
Sub Saharan Africa  |   2.078329   1.326977     1.15   0.252     .5946165    7.264265
     Latin America  |   .6350325   .3074813    -0.94   0.348     .2458363    1.640385
--------------------+----------------------------------------------------------------
              /cut1 |  -4.457378   3.050337                     -10.43593    1.521173
              /cut2 |  -2.838361   3.049647                     -8.815559    3.138837
              /cut3 |  -1.108958   3.049097                     -7.085077    4.867162
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   1.049933   .2774822                      .6254623     1.76247
-------------------------------------------------------------------------------------
Note: Estimates are transformed only in the first equation to odds ratios.
LR test vs. ologit model: chibar2(01) = 1719.82       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class polit
> ics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, or nested dec
> (3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood =  -18951.51  
Iteration 2:  Log likelihood = -18935.781  
Iteration 3:  Log likelihood = -18935.759  
Iteration 4:  Log likelihood = -18935.759  

Refining starting values:

Grid node 0:  Log likelihood = -18161.658

Fitting full model:

Iteration 0:  Log likelihood = -18161.658  
Iteration 1:  Log likelihood =  -18136.07  
Iteration 2:  Log likelihood = -18084.975  
Iteration 3:  Log likelihood = -18062.489  
Iteration 4:  Log likelihood = -18061.514  
Iteration 5:  Log likelihood = -18061.468  
Iteration 6:  Log likelihood = -18061.468  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1631.36
Log likelihood = -18061.468                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Odds ratio   Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
       vfundingx_l1 |   .6364668   .0482676    -5.96   0.000      .548559     .738462
             mxx_l1 |   .9904379   .0084187    -1.13   0.258     .9740743    1.007076
      church_attend |   1.237568   .0093855    28.11   0.000     1.219309      1.2561
          religious |   1.892219   .0835952    14.44   0.000     1.735269    2.063366
            atheist |   .7961393   .1233121    -1.47   0.141     .5876903    1.078523
     satisfied_life |   1.044944   .0075796     6.06   0.000     1.030194    1.059906
              edu08 |    .929855   .0073421    -9.21   0.000     .9155755    .9443572
       social_class |   1.000416    .016402     0.03   0.980     .9687794    1.033085
       politics_l2r |   1.050158   .0070371     7.30   0.000     1.036456    1.064042
          zpcgdplog |   .6832812   .1051495    -2.47   0.013     .5053693    .9238258
           zpolity2 |   .9614902   .0116485    -3.24   0.001     .9389286     .984594
           zdurable |   1.008183   .0048018     1.71   0.087     .9988157    1.017639
          ezmajpctx |   .9886068   .0107885    -1.05   0.294     .9676862     1.00998
            zpoplog |   1.014323   .1283594     0.11   0.911     .7915155     1.29985
    fraction_tc_min |   .0120547   .0138778    -3.84   0.000     .0012625    .1151047
                    |
             region |
       Former USSR  |   1.109935    .598718     0.19   0.847     .3856068    3.194849
              Asia  |   1.259521   1.562188     0.19   0.852     .1107789    14.32035
Sub Saharan Africa  |   1.057796    .714751     0.08   0.934     .2813498    3.977016
     Latin America  |   .7984063   .4142956    -0.43   0.664     .2887572    2.207573
--------------------+----------------------------------------------------------------
              /cut1 |  -8.444036   3.246118                     -14.80631   -2.081762
              /cut2 |  -6.822377   3.245432                     -13.18331   -.4614475
              /cut3 |  -5.090229   3.244721                     -11.44976    1.269308
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|    1.24168   .3238958                      .7446828     2.07037
-------------------------------------------------------------------------------------
Note: Estimates are transformed only in the first equation to odds ratios.
LR test vs. ologit model: chibar2(01) = 1748.58       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. asdoc meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_clas
> s politics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:, or nes
> ted dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) stat(ch2 p ll) append

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood = -18943.101  
Iteration 2:  Log likelihood = -18927.443  
Iteration 3:  Log likelihood =  -18927.42  
Iteration 4:  Log likelihood =  -18927.42  

Refining starting values:

Grid node 0:  Log likelihood = -18159.397

Fitting full model:

Iteration 0:  Log likelihood = -18159.397  
Iteration 1:  Log likelihood = -18147.577  
Iteration 2:  Log likelihood =   -18117.6  
Iteration 3:  Log likelihood = -18087.999  
Iteration 4:  Log likelihood = -18079.613  
Iteration 5:  Log likelihood = -18076.489  
Iteration 6:  Log likelihood =  -18076.48  
Iteration 7:  Log likelihood =  -18076.48  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1614.85
Log likelihood = -18076.48                      Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Odds ratio   Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |   1.316892   .1119097     3.24   0.001     1.114846    1.555555
             mxx_l1 |   .9787905   .0077931    -2.69   0.007     .9636348    .9941846
      church_attend |   1.239602    .009397    28.33   0.000      1.22132    1.258157
          religious |   1.895383   .0837216    14.48   0.000     1.738194    2.066787
            atheist |   .8012484   .1240402    -1.43   0.152     .5915532    1.085277
     satisfied_life |   1.046111    .007584     6.22   0.000     1.031352    1.061081
              edu08 |   .9306107   .0073378    -9.12   0.000     .9163395    .9451043
       social_class |   1.001594   .0164122     0.10   0.923     .9699375    1.034283
       politics_l2r |    1.05115   .0070388     7.45   0.000     1.037444    1.065037
          zpcgdplog |   .8097462    .110696    -1.54   0.123     .6194209    1.058552
           zpolity2 |   .9601811   .0115484    -3.38   0.001     .9378114    .9830844
           zdurable |   .9983659   .0037927    -0.43   0.667     .9909598    1.005827
          ezmajpctx |   1.005899   .0101584     0.58   0.560      .986185    1.026008
            zpoplog |   1.001621    .111577     0.01   0.988     .8051605    1.246018
    fraction_tc_min |   .2033549   .2086943    -1.55   0.121     .0272084    1.519867
                    |
             region |
       Former USSR  |   .7183804   .3559294    -0.67   0.504     .2720304    1.897106
              Asia  |    1.18936   1.297768     0.16   0.874     .1401283    10.09486
Sub Saharan Africa  |   1.169868   .7208949     0.25   0.799     .3496311    3.914384
     Latin America  |   .5361765   .2479761    -1.35   0.178     .2165888    1.327332
--------------------+----------------------------------------------------------------
              /cut1 |  -3.500744   2.794174                     -8.977225    1.975736
              /cut2 |  -1.881682   2.793629                     -7.357095     3.59373
              /cut3 |   -.153785   2.793213                     -5.628382    5.320812
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .9312984   .2321379                      .5713683    1.517964
-------------------------------------------------------------------------------------
Note: Estimates are transformed only in the first equation to odds ratios.
LR test vs. ologit model: chibar2(01) = 1701.88       Prob >= chibar2 = 0.0000
Click to Open File:  Myfile.doc

. 
. *****************************

. * Marginal Effects Analyses *

. *****************************

. 
. meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_
> l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85630.315  
Iteration 2:  Log likelihood = -85372.594  
Iteration 3:  Log likelihood = -85371.867  
Iteration 4:  Log likelihood = -85371.867  

Refining starting values:

Grid node 0:  Log likelihood = -82754.802

Fitting full model:

Iteration 0:  Log likelihood = -82754.802  (not concave)
Iteration 1:  Log likelihood = -82748.023  (not concave)
Iteration 2:  Log likelihood = -82741.652  (backed up)
Iteration 3:  Log likelihood = -82669.758  
Iteration 4:  Log likelihood = -82569.519  
Iteration 5:  Log likelihood = -82469.532  
Iteration 6:  Log likelihood =  -82441.29  
Iteration 7:  Log likelihood = -82430.806  
Iteration 8:  Log likelihood =  -82430.75  
Iteration 9:  Log likelihood =  -82430.75  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   11989.14
Log likelihood = -82430.75                      Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
     lz06funding_l1 |   .0281889   .0139422     2.02   0.043     .0008626    .0555151
             mxx_l1 |   .0058173   .0030766     1.89   0.059    -.0002127    .0118474
      church_attend |    .288619   .0037609    76.74   0.000     .2812479    .2959902
          religious |   .8777018   .0198287    44.26   0.000     .8388383    .9165654
            atheist |  -.6358739   .0699162    -9.09   0.000    -.7729071   -.4988407
     satisfied_life |   .0311928    .003349     9.31   0.000     .0246289    .0377567
              edu08 |  -.0913121   .0035866   -25.46   0.000    -.0983417   -.0842825
       social_class |  -.0132357   .0081307    -1.63   0.104    -.0291715    .0027001
       politics_l2r |   .0471544   .0031741    14.86   0.000     .0409333    .0533754
          zpcgdplog |   .1628553   .0771961     2.11   0.035     .0115538    .3141569
           zpolity2 |   .0031426   .0049527     0.63   0.526    -.0065644    .0128497
           zdurable |  -.0078576   .0026374    -2.98   0.003    -.0130268   -.0026884
          ezmajpctx |  -.0178103   .0046691    -3.81   0.000    -.0269615   -.0086592
            zpoplog |  -.2454276   .0957104    -2.56   0.010    -.4330165   -.0578387
    fraction_tc_min |  -3.083792    .417558    -7.39   0.000    -3.902191   -2.265394
                    |
             region |
       Former USSR  |  -.6844166   .3845175    -1.78   0.075    -1.438057     .069224
              Asia  |   1.391645   .9954373     1.40   0.162     -.559376    3.342666
Sub Saharan Africa  |   .9439547   .4651709     2.03   0.042     .0322366    1.855673
     Latin America  |   .1249988   .3930263     0.32   0.750    -.6453186    .8953163
--------------------+----------------------------------------------------------------
              /cut1 |  -6.940299   1.939921                     -10.74247   -3.138124
              /cut2 |  -4.664661   1.939713                      -8.46643   -.8628929
              /cut3 |  -2.622667   1.939592                     -6.424198    1.178864
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .8665677   .2034077                      .5470181    1.372787
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 5882.23       Prob >= chibar2 = 0.0000

. 
. margins, dydx(lz06funding_l1) post

Average marginal effects                                Number of obs = 78,132
Model VCE: OIM

dy/dx wrt: lz06funding_l1

1._predict: Marginal predicted mean (1.confidence_church), predict(pr outcome(1))
2._predict: Marginal predicted mean (2.confidence_church), predict(pr outcome(2))
3._predict: Marginal predicted mean (3.confidence_church), predict(pr outcome(3))
4._predict: Marginal predicted mean (4.confidence_church), predict(pr outcome(4))

--------------------------------------------------------------------------------
               |            Delta-method
               |      dy/dx   std. err.      z    P>|z|     [95% conf. interval]
---------------+----------------------------------------------------------------
lz06funding_l1 |
      _predict |
            1  |  -.0015509   .0007859    -1.97   0.048    -.0030912   -.0000106
            2  |  -.0026388   .0012978    -2.03   0.042    -.0051826   -.0000951
            3  |  -.0004763   .0003449    -1.38   0.167    -.0011524    .0001997
            4  |   .0046661   .0022998     2.03   0.042     .0001585    .0091737
--------------------------------------------------------------------------------

. 
. est store lz06_maj

. 
. meologit confidence_church lz06funding_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_
> l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood = -18949.878  
Iteration 2:  Log likelihood = -18933.967  
Iteration 3:  Log likelihood = -18933.944  
Iteration 4:  Log likelihood = -18933.944  

Refining starting values:

Grid node 0:  Log likelihood = -18157.202

Fitting full model:

Iteration 0:  Log likelihood = -18157.202  
Iteration 1:  Log likelihood = -18136.124  
Iteration 2:  Log likelihood = -18092.148  
Iteration 3:  Log likelihood = -18074.328  
Iteration 4:  Log likelihood = -18074.038  
Iteration 5:  Log likelihood = -18074.033  
Iteration 6:  Log likelihood = -18074.033  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1616.00
Log likelihood = -18074.033                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
     lz06funding_l1 |  -.1001633   .0258016    -3.88   0.000    -.1507334   -.0495931
             mxx_l1 |  -.0171478   .0081811    -2.10   0.036    -.0331824   -.0011131
      church_attend |   .2136186   .0075845    28.17   0.000     .1987532     .228484
          religious |   .6381599   .0441664    14.45   0.000     .5515953    .7247244
            atheist |  -.2201592   .1548197    -1.42   0.155    -.5236003    .0832818
     satisfied_life |   .0439243   .0072548     6.05   0.000     .0297051    .0581435
              edu08 |  -.0725832     .00789    -9.20   0.000    -.0880473    -.057119
       social_class |   .0010525   .0163898     0.06   0.949    -.0310708    .0331759
       politics_l2r |   .0494329   .0066979     7.38   0.000     .0363054    .0625605
          zpcgdplog |  -.1241533   .1433751    -0.87   0.387    -.4051633    .1568567
           zpolity2 |  -.0377301   .0120816    -3.12   0.002    -.0614096   -.0140506
           zdurable |   .0032531   .0043667     0.74   0.456    -.0053055    .0118117
          ezmajpctx |  -.0020879   .0104645    -0.20   0.842    -.0225979    .0184222
            zpoplog |  -.0238505   .1188376    -0.20   0.841     -.256768    .2090669
    fraction_tc_min |  -2.915481    1.09117    -2.67   0.008    -5.054136   -.7768263
                    |
             region |
       Former USSR  |   .1814228    .516278     0.35   0.725    -.8304635    1.193309
              Asia  |   .7023653    1.15021     0.61   0.541    -1.552004    2.956735
Sub Saharan Africa  |   .7315643   .6384826     1.15   0.252    -.5198385    1.982967
     Latin America  |  -.4540791   .4841978    -0.94   0.348    -1.403089    .4949311
--------------------+----------------------------------------------------------------
              /cut1 |  -4.457378   3.050337                     -10.43593    1.521173
              /cut2 |  -2.838361   3.049647                     -8.815559    3.138837
              /cut3 |  -1.108958   3.049097                     -7.085077    4.867162
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   1.049933   .2774822                      .6254623     1.76247
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 1719.82       Prob >= chibar2 = 0.0000

. 
. margins, dydx(lz06funding_l1) post

Average marginal effects                                Number of obs = 15,879
Model VCE: OIM

dy/dx wrt: lz06funding_l1

1._predict: Marginal predicted mean (1.confidence_church), predict(pr outcome(1))
2._predict: Marginal predicted mean (2.confidence_church), predict(pr outcome(2))
3._predict: Marginal predicted mean (3.confidence_church), predict(pr outcome(3))
4._predict: Marginal predicted mean (4.confidence_church), predict(pr outcome(4))

--------------------------------------------------------------------------------
               |            Delta-method
               |      dy/dx   std. err.      z    P>|z|     [95% conf. interval]
---------------+----------------------------------------------------------------
lz06funding_l1 |
      _predict |
            1  |    .008951   .0023828     3.76   0.000     .0042809    .0136212
            2  |    .007111   .0019632     3.62   0.000     .0032632    .0109588
            3  |   .0006623   .0011827     0.56   0.576    -.0016558    .0029803
            4  |  -.0167243   .0044013    -3.80   0.000    -.0253508   -.0080979
--------------------------------------------------------------------------------

. 
. est store lz06_min

. 
. coefplot (lz06_maj, label(Majority) offset(-0.001)) (lz06_min, label(Minority) offset(-0.001)), vert recast(line) ciopts(rec
> ast(rarea)) graphregion(color(white)) bgcolor(white) yline(0) xlabel( 1 "1" 2 "2" 3 "3" 4 "4") xtitle(Confidence in the Chur
> ch) ytitle(Avg Effects on Marginal Probability) legend(pos(6) row(1)) title(General Religious Funding) saving(lz06)
file lz06.gph saved

. 
. 
. 
. meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2
> r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85578.777  
Iteration 2:  Log likelihood = -85320.969  
Iteration 3:  Log likelihood = -85320.236  
Iteration 4:  Log likelihood = -85320.236  

Refining starting values:

Grid node 0:  Log likelihood = -82751.627

Fitting full model:

Iteration 0:  Log likelihood = -82751.627  (not concave)
Iteration 1:  Log likelihood = -82744.797  (not concave)
Iteration 2:  Log likelihood = -82741.632  (backed up)
Iteration 3:  Log likelihood = -82667.516  
Iteration 4:  Log likelihood = -82565.337  
Iteration 5:  Log likelihood = -82462.116  
Iteration 6:  Log likelihood = -82432.373  
Iteration 7:  Log likelihood = -82421.263  
Iteration 8:  Log likelihood = -82421.188  
Iteration 9:  Log likelihood = -82421.188  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   12005.87
Log likelihood = -82421.188                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
       vfundingx_l1 |  -.2543299   .0540828    -4.70   0.000    -.3603302   -.1483297
             mxx_l1 |   .0068141   .0030977     2.20   0.028     .0007428    .0128854
      church_attend |   .2885299    .003761    76.72   0.000     .2811584    .2959014
          religious |   .8761553   .0198312    44.18   0.000     .8372868    .9150237
            atheist |  -.6417827   .0699076    -9.18   0.000    -.7787991   -.5047664
     satisfied_life |   .0312774   .0033489     9.34   0.000     .0247138    .0378411
              edu08 |  -.0912809   .0035869   -25.45   0.000    -.0983111   -.0842507
       social_class |  -.0138862   .0081324    -1.71   0.088    -.0298254    .0020529
       politics_l2r |   .0470009   .0031742    14.81   0.000     .0407796    .0532223
          zpcgdplog |   .1616952   .0767292     2.11   0.035     .0113087    .3120817
           zpolity2 |   .0036558   .0049453     0.74   0.460    -.0060368    .0133484
           zdurable |  -.0068859   .0026172    -2.63   0.009    -.0120154   -.0017563
          ezmajpctx |  -.0180203   .0045659    -3.95   0.000    -.0269693   -.0090712
            zpoplog |  -.2302189      .0995    -2.31   0.021    -.4252353   -.0352026
    fraction_tc_min |    -3.2902   .4121777    -7.98   0.000    -4.098054   -2.482347
                    |
             region |
       Former USSR  |  -.4243227    .392082    -1.08   0.279    -1.192789     .344144
              Asia  |   1.463924   1.040657     1.41   0.160    -.5757258    3.503574
Sub Saharan Africa  |   .6737275   .4890047     1.38   0.168    -.2847042    1.632159
     Latin America  |   .3066156   .4098282     0.75   0.454     -.496633    1.109864
--------------------+----------------------------------------------------------------
              /cut1 |  -7.533581   2.001203                     -11.45587   -3.611296
              /cut2 |  -5.257719   2.001004                     -9.179615   -1.335823
              /cut3 |  -3.215087   2.000873                     -7.136727    .7065522
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .9528289   .2207482                      .6050781    1.500439
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 5798.10       Prob >= chibar2 = 0.0000

. 
. margins, dydx(vfundingx_l1) post

Average marginal effects                                Number of obs = 78,132
Model VCE: OIM

dy/dx wrt: vfundingx_l1

1._predict: Marginal predicted mean (1.confidence_church), predict(pr outcome(1))
2._predict: Marginal predicted mean (2.confidence_church), predict(pr outcome(2))
3._predict: Marginal predicted mean (3.confidence_church), predict(pr outcome(3))
4._predict: Marginal predicted mean (4.confidence_church), predict(pr outcome(4))

------------------------------------------------------------------------------
             |            Delta-method
             |      dy/dx   std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
vfundingx_l1 |
    _predict |
          1  |   .0141301   .0034236     4.13   0.000     .0074199    .0208403
          2  |   .0236188   .0049371     4.78   0.000     .0139423    .0332953
          3  |   .0038864   .0024326     1.60   0.110    -.0008813    .0086542
          4  |  -.0416353   .0086343    -4.82   0.000    -.0585583   -.0247124
------------------------------------------------------------------------------

. 
. est store vfund_maj

. 
. meologit confidence_church vfundingx_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class politics_l2
> r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood =  -18951.51  
Iteration 2:  Log likelihood = -18935.781  
Iteration 3:  Log likelihood = -18935.759  
Iteration 4:  Log likelihood = -18935.759  

Refining starting values:

Grid node 0:  Log likelihood = -18161.658

Fitting full model:

Iteration 0:  Log likelihood = -18161.658  
Iteration 1:  Log likelihood =  -18136.07  
Iteration 2:  Log likelihood = -18084.975  
Iteration 3:  Log likelihood = -18062.489  
Iteration 4:  Log likelihood = -18061.514  
Iteration 5:  Log likelihood = -18061.468  
Iteration 6:  Log likelihood = -18061.468  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1631.36
Log likelihood = -18061.468                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
       vfundingx_l1 |  -.4518231   .0758368    -5.96   0.000    -.6004605   -.3031857
             mxx_l1 |  -.0096081   .0084999    -1.13   0.258    -.0262677    .0070514
      church_attend |    .213148   .0075838    28.11   0.000      .198284     .228012
          religious |   .6377505   .0441784    14.44   0.000     .5511624    .7243385
            atheist |  -.2279811   .1548875    -1.47   0.141    -.5315551    .0755929
     satisfied_life |   .0439635   .0072536     6.06   0.000     .0297468    .0581802
              edu08 |  -.0727266    .007896    -9.21   0.000    -.0882024   -.0572508
       social_class |   .0004156   .0163952     0.03   0.980    -.0317184    .0325495
       politics_l2r |   .0489411    .006701     7.30   0.000     .0358073    .0620748
          zpcgdplog |  -.3808488   .1538891    -2.47   0.013    -.6824658   -.0792318
           zpolity2 |  -.0392709    .012115    -3.24   0.001    -.0630159   -.0155259
           zdurable |   .0081499   .0047628     1.71   0.087     -.001185    .0174849
          ezmajpctx |  -.0114586   .0109129    -1.05   0.294    -.0328474    .0099302
            zpoplog |   .0142215   .1265469     0.11   0.911    -.2338059    .2622488
    fraction_tc_min |  -4.418305   1.151241    -3.84   0.000    -6.674696   -2.161913
                    |
             region |
       Former USSR  |   .1043013   .5394173     0.19   0.847    -.9529371     1.16154
              Asia  |   .2307314   1.240303     0.19   0.852    -2.200219    2.661682
Sub Saharan Africa  |   .0561877   .6756982     0.08   0.934    -1.268156    1.380532
     Latin America  |  -.2251377   .5189032    -0.43   0.664    -1.242169    .7918939
--------------------+----------------------------------------------------------------
              /cut1 |  -8.444036   3.246118                     -14.80631   -2.081762
              /cut2 |  -6.822377   3.245432                     -13.18331   -.4614475
              /cut3 |  -5.090229   3.244721                     -11.44976    1.269308
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|    1.24168   .3238958                      .7446828     2.07037
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 1748.58       Prob >= chibar2 = 0.0000

. 
. margins, dydx(vfundingx_l1) post

Average marginal effects                                Number of obs = 15,879
Model VCE: OIM

dy/dx wrt: vfundingx_l1

1._predict: Marginal predicted mean (1.confidence_church), predict(pr outcome(1))
2._predict: Marginal predicted mean (2.confidence_church), predict(pr outcome(2))
3._predict: Marginal predicted mean (3.confidence_church), predict(pr outcome(3))
4._predict: Marginal predicted mean (4.confidence_church), predict(pr outcome(4))

------------------------------------------------------------------------------
             |            Delta-method
             |      dy/dx   std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
vfundingx_l1 |
    _predict |
          1  |   .0367543   .0069076     5.32   0.000     .0232157    .0502929
          2  |   .0312497      .0047     6.65   0.000     .0220378    .0404615
          3  |   .0072151     .00541     1.33   0.182    -.0033884    .0178185
          4  |   -.075219   .0118656    -6.34   0.000    -.0984751    -.051963
------------------------------------------------------------------------------

. 
. est store vfund_min

. 
. coefplot (vfund_maj, label(Majority) offset(-0.001)) (vfund_min, label(Minority) offset(-0.001)), vert recast(line) ciopts(r
> ecast(rarea)) graphregion(color(white)) bgcolor(white) yline(0) xlabel( 1 "1" 2 "2" 3 "3" 4 "4") xtitle(Confidence in the Ch
> urch) ytitle(Avg Effects on Marginal Probability) legend(pos(6) row(1)) title(Exclusive Religious Funding) saving(vfund)
file vfund.gph saved

. 
. 
. 
. meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class poli
> tics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if majority_fixed==1 || ccode:

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -95590.536  
Iteration 1:  Log likelihood = -85642.458  
Iteration 2:  Log likelihood = -85385.817  
Iteration 3:  Log likelihood = -85385.081  
Iteration 4:  Log likelihood = -85385.081  

Refining starting values:

Grid node 0:  Log likelihood =  -82748.28

Fitting full model:

Iteration 0:  Log likelihood =  -82748.28  (not concave)
Iteration 1:  Log likelihood = -82741.473  (backed up)
Iteration 2:  Log likelihood = -82736.726  (backed up)
Iteration 3:  Log likelihood = -82664.675  
Iteration 4:  Log likelihood = -82565.153  
Iteration 5:  Log likelihood = -82465.497  
Iteration 6:  Log likelihood = -82437.267  
Iteration 7:  Log likelihood = -82426.798  
Iteration 8:  Log likelihood = -82426.741  
Iteration 9:  Log likelihood = -82426.741  

Mixed-effects ologit regression                 Number of obs     =     78,132
Group variable: ccode                           Number of groups  =         51

                                                Obs per group:
                                                              min =        174
                                                              avg =    1,532.0
                                                              max =      6,526

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =   11995.30
Log likelihood = -82426.741                     Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |  -.0927677   .0266623    -3.48   0.001    -.1450249   -.0405104
             mxx_l1 |   .0045935   .0030939     1.48   0.138    -.0014705    .0106575
      church_attend |    .288736    .003761    76.77   0.000     .2813646    .2961073
          religious |   .8769958   .0198294    44.23   0.000     .8381308    .9158607
            atheist |   -.645832   .0699536    -9.23   0.000    -.7829386   -.5087254
     satisfied_life |   .0311186   .0033487     9.29   0.000     .0245552     .037682
              edu08 |  -.0915676   .0035877   -25.52   0.000    -.0985994   -.0845359
       social_class |  -.0137969   .0081319    -1.70   0.090     -.029735    .0021413
       politics_l2r |   .0470076   .0031743    14.81   0.000     .0407861    .0532291
          zpcgdplog |   .2238082   .0761193     2.94   0.003     .0746172    .3729993
           zpolity2 |   .0024269   .0049543     0.49   0.624    -.0072833     .012137
           zdurable |  -.0072119   .0025872    -2.79   0.005    -.0122827   -.0021412
          ezmajpctx |  -.0157509   .0046241    -3.41   0.001     -.024814   -.0066878
            zpoplog |  -.2814766   .0985131    -2.86   0.004    -.4745587   -.0883946
    fraction_tc_min |  -3.194759   .4174929    -7.65   0.000     -4.01303   -2.376488
                    |
             region |
       Former USSR  |  -.4581258    .381864    -1.20   0.230    -1.206566     .290314
              Asia  |   1.649603   1.002732     1.65   0.100    -.3157154    3.614922
Sub Saharan Africa  |   1.207318   .4690428     2.57   0.010      .288011    2.126625
     Latin America  |   .2129496   .3949811     0.54   0.590    -.5611991    .9870983
--------------------+----------------------------------------------------------------
              /cut1 |  -6.939627   1.981942                     -10.82416   -3.055092
              /cut2 |  -4.663633   1.981733                     -8.547758   -.7795074
              /cut3 |  -2.621301   1.981599                     -6.505164    1.262562
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .8819603    .207134                      .5565945    1.397524
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 5916.68       Prob >= chibar2 = 0.0000

. 
. margins, dydx(lz07entanglement_l1) post

Average marginal effects                                Number of obs = 78,132
Model VCE: OIM

dy/dx wrt: lz07entanglement_l1

1._predict: Marginal predicted mean (1.confidence_church), predict(pr outcome(1))
2._predict: Marginal predicted mean (2.confidence_church), predict(pr outcome(2))
3._predict: Marginal predicted mean (3.confidence_church), predict(pr outcome(3))
4._predict: Marginal predicted mean (4.confidence_church), predict(pr outcome(4))

-------------------------------------------------------------------------------------
                    |            Delta-method
                    |      dy/dx   std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |
           _predict |
                 1  |   .0051369   .0015703     3.27   0.001     .0020592    .0082146
                 2  |   .0085781   .0024562     3.49   0.000     .0037641    .0133922
                 3  |     .00158    .000943     1.68   0.094    -.0002683    .0034283
                 4  |   -.015295   .0043845    -3.49   0.000    -.0238884   -.0067016
-------------------------------------------------------------------------------------

. 
. est store lz07_maj

. 
. meologit confidence_church lz07entanglement_l1 mxx_l1 church_attend religious atheist satisfied_life edu08 social_class poli
> tics_l2r zpcgdplog zpolity2 zdurable ezmajpctx zpoplog fraction_tc_min i.region if minority_fixed==1 || ccode:

Fitting fixed-effects model:

Iteration 0:  Log likelihood = -20200.116  
Iteration 1:  Log likelihood = -18943.101  
Iteration 2:  Log likelihood = -18927.443  
Iteration 3:  Log likelihood =  -18927.42  
Iteration 4:  Log likelihood =  -18927.42  

Refining starting values:

Grid node 0:  Log likelihood = -18159.397

Fitting full model:

Iteration 0:  Log likelihood = -18159.397  
Iteration 1:  Log likelihood = -18147.577  
Iteration 2:  Log likelihood =   -18117.6  
Iteration 3:  Log likelihood = -18087.999  
Iteration 4:  Log likelihood = -18079.613  
Iteration 5:  Log likelihood = -18076.489  
Iteration 6:  Log likelihood =  -18076.48  
Iteration 7:  Log likelihood =  -18076.48  

Mixed-effects ologit regression                 Number of obs     =     15,879
Group variable: ccode                           Number of groups  =         50

                                                Obs per group:
                                                              min =          4
                                                              avg =      317.6
                                                              max =      2,386

Integration method: mvaghermite                 Integration pts.  =          7

                                                Wald chi2(19)     =    1614.85
Log likelihood = -18076.48                      Prob > chi2       =     0.0000
-------------------------------------------------------------------------------------
  confidence_church | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |   .2752741   .0849802     3.24   0.001     .1087159    .4418323
             mxx_l1 |  -.0214377    .007962    -2.69   0.007    -.0370429   -.0058324
      church_attend |   .2147902   .0075806    28.33   0.000     .1999324     .229648
          religious |   .6394209   .0441713    14.48   0.000     .5528466    .7259951
            atheist |  -.2215843   .1548086    -1.43   0.152    -.5250036    .0818351
     satisfied_life |   .0450795   .0072497     6.22   0.000     .0308704    .0592886
              edu08 |  -.0719142   .0078849    -9.12   0.000    -.0873684     -.05646
       social_class |   .0015925   .0163861     0.10   0.923    -.0305237    .0337086
       politics_l2r |   .0498847   .0066963     7.45   0.000     .0367603    .0630091
          zpcgdplog |  -.2110344   .1367045    -1.54   0.123    -.4789703    .0569016
           zpolity2 |  -.0406334   .0120273    -3.38   0.001    -.0642065   -.0170603
           zdurable |  -.0016355    .003799    -0.43   0.667    -.0090813    .0058103
          ezmajpctx |   .0058819   .0100988     0.58   0.560    -.0139113    .0256752
            zpoplog |   .0016194   .1113965     0.01   0.988    -.2167137    .2199525
    fraction_tc_min |  -1.592802   1.026256    -1.55   0.121    -3.604228    .4186231
                    |
             region |
       Former USSR  |   -.330756   .4954609    -0.67   0.504    -1.301841    .6403294
              Asia  |   .1734151   1.091148     0.16   0.874    -1.965197    2.312027
Sub Saharan Africa  |   .1568907   .6162192     0.25   0.799    -1.050877    1.364658
     Latin America  |  -.6232918   .4624896    -1.35   0.178    -1.529755    .2831711
--------------------+----------------------------------------------------------------
              /cut1 |  -3.500744   2.794174                     -8.977225    1.975736
              /cut2 |  -1.881682   2.793629                     -7.357095     3.59373
              /cut3 |   -.153785   2.793213                     -5.628382    5.320812
--------------------+----------------------------------------------------------------
ccode               |
          var(_cons)|   .9312984   .2321379                      .5713683    1.517964
-------------------------------------------------------------------------------------
LR test vs. ologit model: chibar2(01) = 1701.88       Prob >= chibar2 = 0.0000

. 
. margins, dydx(lz07entanglement_l1) post

Average marginal effects                                Number of obs = 15,879
Model VCE: OIM

dy/dx wrt: lz07entanglement_l1

1._predict: Marginal predicted mean (1.confidence_church), predict(pr outcome(1))
2._predict: Marginal predicted mean (2.confidence_church), predict(pr outcome(2))
3._predict: Marginal predicted mean (3.confidence_church), predict(pr outcome(3))
4._predict: Marginal predicted mean (4.confidence_church), predict(pr outcome(4))

-------------------------------------------------------------------------------------
                    |            Delta-method
                    |      dy/dx   std. err.      z    P>|z|     [95% conf. interval]
--------------------+----------------------------------------------------------------
lz07entanglement_l1 |
           _predict |
                 1  |  -.0259118   .0084068    -3.08   0.002    -.0423889   -.0094347
                 2  |  -.0193148   .0064219    -3.01   0.003    -.0319015   -.0067281
                 3  |  -.0002058   .0028488    -0.07   0.942    -.0057894    .0053777
                 4  |   .0454324   .0141788     3.20   0.001     .0176424    .0732224
-------------------------------------------------------------------------------------

. 
. est store lz07_min

. 
. coefplot (lz07_maj, label(Majority) offset(-0.001)) (lz07_min, label(Minority) offset(-0.001)), vert recast(line) ciopts(rec
> ast(rarea)) graphregion(color(white)) bgcolor(white) yline(0) xlabel( 1 "1" 2 "2" 3 "3" 4 "4") xtitle(Confidence in the Chur
> ch) ytitle(Avg Effects on Marginal Probability) legend(pos(6) row(1)) title(Religious Entanglement) saving(lz07)
file lz07.gph saved

. 
. gr combine lz06.gph vfund.gph lz07.gph, ycommon row(1) ysize(2) xsize(6) scale(1)

. log close
      name:  <unnamed>
       log:  /dataverse files/Re-Assessing Governmental Corruption of Religion-2026-logfile.log
  log type:  text
 closed on:  29 Jan 2026, 10:21:18
------------------------------------------------------------------------------------------------------------------------------
